NPS-MA-93-007 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


POWER  ITERATIONS  AND  THE  DOMINANT 
EIGENVALUE  PROBLEM 


by 

Jeffery  J.  Leader 
ft 

Technical  Report  For  Period 
March  1992  -  June  1992 


Approved  for  public  release;  distribution  unlimited 

PedDocs       pared  for:   Naval  Postgraduate  School 
D  208.14/2  Monterey,  CA  93943 

NPS-MA-93-007 


£0*.  u\\x  OPS-  mA  -1 V  C»7 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 

Rear  Admiral  T.A.  Mercer  .- .     _.  .,  , 

.  .  Harrison  Shull 

Superintendent  Provost 

This  report  was  prepared  in  conjunction  with  research  conducted 
for  the  Naval  Postgraduate  School  and  funded  by  the  Naval 
Postgraduate  School. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 
This  report  was  prepared  by: 


Jnclassif  ied 


£C..nH  Y  CLASS.F'CATIQN   QF    'H.S 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  93943-5101 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-011 


a   REPORT  SECURITY  CLASSIFICATION 

Jnclassif  ied 


lb    RESTRICTIVE  MARKINGS 


a    SECURITY  CLASSIFICATION  AUTHORITY 


b   DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3     DlSTRIBUTlON/AVAILABlL.TY  OF  REPORT 

Approved  for  public  release 
distribution  unlimited 


.   PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

IPS-MA-93-007 


5    MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

NPS-MA-93-007 


a    NAME  OF  PERFORMING  ORGANIZATION 

laval   Postgraduate   School 


6b    OFFICE  SYMBOL 
(If  applicable) 

MA 


7a.  NAME  OF  MONITORING  ORGANIZATION 


Naval   Postgraduate   School 


c.  ADDRESS  [City,  State,  and  ZIP  Code) 
'aval    Postgraduate    School 

[onterey,    CA  93943 


7b  ADDRESS  (C/'ry.  State,  and  ZiPCode) 
Naval  Postgraduate  School 
Monterey,    CA   93943 


a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

aval  Postgraduate  School 


8b    OFFICE  SYMBOL 
(If  applicable) 

MA 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


O&MN  Direct  Funding 


c.  ADDRESS  (City,  State,  and  ZIP  Code) 

aval  Postgraduate  School 
onterey,  CA  93943 


:0    SOURCE  OF  FUNDING  NUMBE'S 


PROGRAM 
ELEMENT  NO 


PROJECT 
NO 


TASK 
NO 


WORK  UNIT 
ACCESSION  NO 


1.  TITLE  (Include  Security  Classification) 

ower  Iterations  and  the  Dominant  Eigenvalue  Problem 


2.  PERSONAL  AUTHOR(S) 

effery   J.    Leader 


3a    TYPE  OF  REPORT 

echnical   Report 


3b    TIME  COVERED 

from  3-92 to  6-92 


14    DATE  OF  REPORT    {Year,  Month.  Day) 

12-15-92 


5    PAGE  COU.V 


11 


5   SUPPLEMENTARY  NOTATION 


COSATI  CODES 


FIELD 


GROUP 


SUB-GROUP 


18    SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  numoer) 

Power  iteration,  Power  method 


J   ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  orbits  of  an  iterative  numerical  method  for  the  dominant  eigenvalue  problem  are 
lalyzed  from  a  discrete  dynamical  systems  perspective.   It  is  shown  that  the  method 
m  extract  more  information  thant  the  standard  power  method  but  at  greater  computational 

)St. 


D    DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
LjJuNCLASSlFlED/UNLIMITED       □  SAME  AS  RPT  □  DTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


la    NAME  OF  RESPONSIBLE   INDIVIDUAL 

Jeffery  J.  Leader 


22b  TELEPHONE  (Include  Area  Code) 


Ut0m%>=22£ 


22c    OFFICE  SYMBOL 

MA/Le 


)Form  1473,   JUN  86 


Previous  editions  are  obsolete 
S/N   0102-LF-014-6603 


SECURITY  CLASSIFICATION  OF   THIS  PAGE 


POWER  ITERATIONS  AND  THE  DOMINANT  EIGENVALUE  PROBLEM 

JEFFERY  J.  LEADEF: 

r^,i —  -  4-  «-  -  <r-  +  t-  +-   m-  .■*  v- -„-_  4  ;  .—  r-    w  2  Va  1   r-.-.»-4.  ^-»^  -  .-i,  T-  +  — .  •—  -i .,— ■ 


>i  i^-x  «—  -^ 


Abstract.   Tht 

dominant  eigenvalue  problem  are  analyzed  from  a  ciscrete 
dynamical  systems  perspective.  It  is  shown  that  the  method 
can  extract  more  information  than  the  standard  power  method 
but  at  greater  computational  cost. 
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to    the     pouter     mj&th.od.    for     finding    the     dominant     eigenvalue     of     a 
real     ma J  rix, 
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iy     slower       than      the      power       method 
can      provide      extra      information      about      the      dominant 


eiaenvaj  ubCs'J       of       a       matrix      in       certain       cases 
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CattractorsD     in    various    cases    Cin    the    spirit    of     [93D. 

Thie    iteration    CI.  ID     is    considered    in    a     different-     context     a: 
a    special     case    of    the    [R   -»[R      map 


V  =    A*V      +    B*V  /IIV    I 

n+l  r>  n  n 


CI  .  3D 

in  [101-  based  on  work  in  [2]  Calso  reported  in  C31D.  Further 
details  on  the  iteration  C 1  .  3D  may  be  found  in  [10,12,133.  and 
t  he  f  or  thcomi ng  [41.  Al though  C 1 . 3D  is  onl y  a  1 i  near 
perturbation  of  the  well -behaved  iteration  CI. ID,  it  exhibits 
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course,  CI.  ID  can  also  be  viewed  as  CI.  S3  with  a  change  of 
nor  mal  i  r.ati  on ,  and  much  is  known  about  the  numerical  method 
given  by  C I . 2D  Csee  also  IS,     p.  3623 D. 

2.  The  Pover  Method 

The  power   method  CI.  2D   has  the  property  that   if   A  is   a 

nondef ecti ve  matrix  with  a  dominant  ei  aer.  value,  say  ?,  ,  and  v 

i  o 

has  a  nonzero  projection  on  an  eigenvector  associated  with 
this    dominant,    eigenvalue,     then 


fj   — >a  dS  n — >oo 
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and  v     converges  to  an  eigenvector  associated  with  X  and  with 

r>  l 

unit  /   norm.    If  u   does  not  have  a  nonzero  component  along 
a  o  f  v 

an  gi  c^'uVcC  ucr  a  s  s  oc  i  a  ^  e*~ <  wi  un  A-  a  no  l  ni  l  ni  t.@  Precision 
cri  thine*  ic  is  used  then  we  must  consider  the  eigenvalue  of 
laraest  modulus  along  which  v>  does  have  a  nonzero  component. 
In  actual  computations,  however,  a  component  along  an 
eigenvector  associated  with  X  would  almost  certainly  be 
introduced  eventually  and  magnified  in  successive  iterations 
CI,  p.  1453.  The  same  results  are  found  if  A  is  defective 
Cconsidering  now  principal  vectors  [7,  p.  3]  rather  than  just 
eigenvectors}  but  the  convergence  is  much  slower. 

When  the  dominant  eigenvalue  is  real,  the  method  converges 
to  a  fixed  point.  However,  when  the  dominant  eigenvalue  is  a 
complex  conjugate  pair,  the  method  generally  fails  to 
converge.  Methods  exist  to  recover  information  in  such  cases 
[6,  p. £57]  but  they  tend  to  be  somewhat  involved. 


Vr'-  -  --       /       -  <-  — ■      i<  i»-c;      »-  — ■-  1  ■  »«-•      +  •—      -       r »-.  -.  t-,  -. .-      .-.♦"      »-  -I    —  ••      ,...-      i —  ,,- 

O  w 

—      A*  v     -    i    v      li  ,  c>  =     A**y  ,-iiA**y      II 

n+i  n  n     00  n  +  1  "  n  +  1  '   n     00 

and  in  this  formulation  v      need  not  be  calculated  until  it  is 

r> 

actually    needed    Cto    estimate    X.  J .        Then    the    iteration    for     v 

1  "n+i 

is    the    siine    as    CI.  II)     except    for     the     particular     /      norm    used    in 

r 

the    normalization.         For     this    reason    we    sometimes     refer     to    the 
quantity    /j      in    C 1  .  2 J     as    the    sif,ed    /      norm. 
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for     every    n>l     and    for     any    V       which    is     nonzero.         Clearly    G    is 

o 

positive  definite  symmetric,  and  so  the  points  V  ,  V  ,  V  ,.  .  . 
must  all  lie  on  the  hyper ell ipse  defined  by 

VT*G*V  =  1  C  3 .  1 5 

in  K  . 

If  B  is  singular,  a  similar  result  holds.  In  order  to 
handle  simultaneously  both  the  case  where  B  is  simple  and  the 
case  where  B  is  defective  we  state  the  result  in  terms  of 
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Theorem  1:  Suppose  Bis  a  real  square  mxm  matrix  with  0<q<rr, 
null  eicenvalues.  If  q> 0  and  V  has  a  nonzero  component  alona 
a  principal  vector  associated  with  a  nonzero  eigenvalue,  then 
all  orb  ts  of  the  power  iteration  CI. 13  are  constrained  to  a 
hyper el  lips©  in  Cm— q3 —dimensions  Cf  or  all  but  finitely  many 
n3  .  Otf  r-rwise,  the  orbit  reaches  the  origin  in  finitely  many 
i  terati ens  . 
Proof:        First ,     note    that 

V      =    B*V   /WW    II 
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be      the      Jordan      normal       form      of       B      for       some      nonsingular       R 

Substituting    this    into    C3.  £3    gives 

V      =    CR*J*P"13r'^V    /IIBr'-1*V    II 
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where    Z   =R     *V      and 
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is     a     scalar     Cfor     each     n3 .  CI  earl  v,      if     Z       has     no     nonzero 
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component     along     a     principal     vector     of     J     C equi val entl  y ,      if     V 

has    no    nonzero    component    along    a    principal     vector     of     B3     that 

is     associated     with     a     non-null     eigenvalue,     the     term     Jr'*Z       in 
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undtf  iiH  d        &.r»d        the        iteration        stops.  Othwr  wise,         for         n 

sufficiently  large  Cit  suffices  that  n>m)  ,  all  Jordan  blocks 
in  J  associated  with  a  null  eigenvalue  will  have  become  blocks 
of  entirely  zeros  in  J  Csince  these  Jordan  blocks  are 
nil  potent}.  Now,  the  principal  vectors  belonging  to  a  given 
Jordan      "lock      do     not      interact      with      the      remaining      principal 

vectors.       in      the      sense      that      if       x>       is      a      principal       vector 
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associated       with       J  CW,        a       Jordan       block       of       J,        then       J  '*u 

v  J 

involves  only  a  linear  combination  of  principal  vectors  of  J 
that,     an       also     associated     with     J  C  \~>  .  Therefore     for     n     large 

enough  that  all  nil  potent  Jordan  blocks  have  become  entirely 
zei  o    submatnces,     the    vector 

jn*z 

o 

can     b^     written     in     terms     of     a     basis     consisting     of      only     the 

remaining  C m-qD  principal  vectors.  Thus  the  iteration  lies  in 
a  C  m-q)  -di  mensi  onal  subspace  of  CP  ,  and  a  suitable  change  of 
variables  c*r  then  be  used  to  transform  the  iteration  into  one 
of     the    :  or  m 

w         =    C*V  x-HW   II 
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where      V        is      a      C  m-qj  -vector      for      every     n      and      C      is      a.      real 

n 

nonsingular  Cm-q) xC m-qD  matrix.  Hence,  in  this  subspace,  the 
iteration  is  constrained  to  the  hyperellipse  determined  by  the 
ma  t  r  i  x 

_T  -1 

CC*C  3 

Cas  was  shown  previously  for  the  nonsingular  case)  for  all  but 
finitely  many  n.  a 

We  emphasize  that  this  is  not  an  asymptotic  result;  after  a 
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hyperell ipse     C assuming      infinite     precision^.  We*      now     wish      to 

look     at    the    orbits    on    the    attract] ng    hyper  el 1 i pses . 

4.     Limit    Orbits    for    Nondefective    Matrices 


Suppose     that      B     is      nondefective     and      nonsingular,       and      lei 
be       a       set       of       ei aenvector s       associated       with       the 
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rhen    any    V   «ef£        may    be    written    as 
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where  a      is  the  component  of  V   along  x  .    Computing  the  power 

V  o  l 

iteration    CI.  ID    by   C3.2D    gives 
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i  cjt     ri-ij.  . 

Now     suppose     that      B     has     a     CrepeatedJ     dominant      figenva.Iu^     X 
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with    multiplicity    lir^m,     that    is. 


and 


A        =     A. 
1  2 
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=      A. 


>     IX 


Removi  no  a  factor  of  X   from  the  numerator  and  X 
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from     the 


denominator    of    C4.13     give; 
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n  1^11 
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•     +    a  x      +    OCX       s\ 
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and    c  1  e.-  r  1  y  ,     i  n    t  he    1  i  mi  t 

V   — >y   x 
where    x    is    a    uni  t    vector    in    span<x   ,  .  .  .  ,  x  >  ,     and 
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nonzero : 


C  pr  ovi  de  j        that        at        least        one        of 

otherwise.      we      beqi n      the      analysis      anew     bv     consider!  nc 
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Hence  i  3  X  >0  the  iteration  tends  to  a  fixed  point,  and  if 
X  <0  the  iteration  tends  to  a  symmetric  Cin  the  origin!)  period 
two  cycle  on  the  points  ±X  x. 

Now  suppose  that  B  has  a  real   positive  eigenvalue  X   of 
mui  1 1  pi i  z 1 1 y  p ,  so  that 
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real     negative    eigenvalue    -X      of     multiplicity    q, 


that 
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p+1        p+2 

and  suppose  further  that 
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Proceeding    as    before,     we    have    that 
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i  =  p  +  l 

,7    given    in    C4.  EjO     where    the    ±    is    positive    when    n    is 


even      and      negative      when      n      is       odd      in       the      numerator,       and 
contrarily     in     the     denominator.  In     the     limit,      we     have     that 
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Since    A      is    positive,     >'    =X      lor     all     ri ,     and    so    we    have    ct    period 
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which    in    general     is    not    a    symmetric    orbit    Cin    the    origin}.     We 

will     not    in    general     havs     II V    II  =X      in    this    case,     but    note    that 
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so    that 
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in  the  limit  of  large  n;  the  modulus  of  the  dominant 
eigenvalue  is  the  geometric  mean  of  the  norms  of  two 
successive  iterates  Cin  the  limits.  Note  that  the  asymmetry 
of  the  orbit  allows  us  to  distinguish  this  case  from  the  case 
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Now  let  us  suppose  that  B  has  a  complex  conjugate  pair  of 
eigenvalues  that  is  dominant,  i.e.   that  X      and  X   are  complex 
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conjugates,     and 
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From    C 4 . 13    we    have 
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where    ^-argCX  J.        In    the    limit,     this    becomes 
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Thus    if     6    is     such    that    expCi0!>=expCiCn+lD03     for     some    n     Ci.e. 

if     X        and     X        are     both     in     [R        for      some     nD      then      we      have      an 
i  i 
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the  orbit  is  aperiodic  on  the  underlying  hyper el 1 i pse.   If  the 

orhjt  is  period  n  then  we  note  that  the  geometric  mean  of  the 

norms  o:   n  consecutive  iterates  tends  to   |  A  |  ,   as  is  easily 

seen  by  writing  out  tne  product  of  n  iterates  and  noting  the 

cancel  1 ati on . 

In     a     similar     way,      if     B     has     two     p^airs     of     complex     conjugate 

eigenvalues     of     equal      modulus     and     all      other      eigenvalues     of     B 

have     lesser      moduli  ,      then     we     have     a     situation     like     that     of 

C4.3D    save    that    there    is    an    additional    angle    to    be    considered. 

Henc<      if     the     firsi      pair     alone     would     give    a     periodic     orbit     of 

period     •,,      and     the     second     pair      alone     would     give     a     periodic 

orbit     of     period     q,      then     the     orbit     of     the     iteration     will      be 

periodic     with     period     n=lcmCp,qj.      If     either     pair      alone     would 

give    an    aperiodic     orbit  ,     then     the    orbit     is     aperiodic.         Again 

the     cieoi  e'.ric     mean     cf      n     consecutive     iterates      tends      to      |  X. 

i 

when      th~      orbit      is      asymptotically      period      n.  The      obvious 

general i lati on  holds  for  more  than  two  complex  conjugate  pairs 
of  equal  moduli.  In  particular,  in  an  even— dimensional  space, 
say  of  dimension  £r  ,  choosing  B  to  be  a  £rx2r  matrix  with  r 
complex  conjugate  pairs  of  eigenvalues,  all  of  equal  modulus 
and    such    that 

\r'    «?   0? 
for       all        i =1 , .  .  .  , £r       and       n=l ,2,3,...        gi  ves       a       method       for 
generating    a    sequence    of     points    on    the    ellipse    given    by    C3.1D 
by    iterating    CI. ID    with    an    arbitrary    nonzero    V    . 

If    in    addition    to    some    number    of    complex    conjugate    pairs    of 
eigenvalues      with      equal       modulus       there      are      some      number       of 


1  n 
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then  it  is  immediate  thai  the  iteration  is  aperiodic  if  any 
compl  e>:  conjugate  eigenvalue  gives  an  aperiodic  orbit,  and 
periodic  otherwise,  with  period  given  by  the  least  common 
multiple  of  the  individual  complex  conjugate  pairs  and  the 
period    tv,'c>    due    to    the    negative    real     eigenvalues,     if    present. 

We    have    established    the    following    theorem: 
Theorfm    L:        Consider     the    iteration 

v       =  b*v  ,■  II  v  i; 

n+l  n  n 

where  &  is  a  real   mxm  nondefective  matrix  and  V   is  a  oi  ven 

o 

nonzero      m-vector  .       Suppose-      that.      V        has      a      nonzero      component 

o 

alone  an  eigenvector  of  &  which  is  associated  with  an 
eigenvalue  of  B  with  maximum  modulus.  Then  all  iterates  for 
n>rn  are  constrained  to  an  ellipse  in  a  subspace  of  CF  "  of 
dimension  rankCBJ,  and  the  asymptotic  behaviour  oi  the 
iterates    is    as     follows: 

i  J     if     B    has     a     multiple     r^^l      dominant     eigenvalue,      then     the 
iteration    tends    to    a    fixed    point; 

ii~>     if     B    has    a    multiple    negative    real     eigenvalue,     then    the 
iteration    tends    to    a    symmetric    period    two    cycle; 

iii'J      if      B      has      both      positive      and      negative      real       dominant 
eigenvalues,     then    the    iteration    tends    to    a    period    two    cycle; 

iv)  if  B  has  a  complex  conjugate  dominant  eigenvalue,  then 
the  iteration  tends  to  a  period  n  cycle  if  the  eigenvalues  are 
a  multiple  of  an  nth  root  of  unity  and  is  aperiodic  otherwise; 
v3>  if  B  has  multiple  complex  conjugate  dominant  eigenvalues 
and  multiple  real  dominant  eigenvalues,  then  the  iteration 
tends      to     a      periodic      orbit      with     period     equal      to      the      least 


11 


ei  genval ues  when  all  eigenvalues  give  rise  to  periodic  orbits, 
and    is    aperiodic    otherwise. 

5.     Limit    Orbits    for    Defective    Matrices 

The       c oove       results       are       essentially       unchanged       if        B       is 

defective.  Let     B     be     defective     and     let     x     be     an     eigenvector 

associated     with     a     dominant     eigenvalue     X.         Let     y     be     a     linear 

combination    of     principal     vectors    of    B    associated    with    X.        Then 

the    i ter ati on 

n 

B  y 

converges  to  x  as  OQlSrO  CIO,  p.  883  ,  [15,  p.  5823  .  Since  the 
principal  vectors  corresponding  to  distinct  eigenvalues  are 
noni  interacting,  it  is  clear  that  the  qualitative  results  of 
Theorem  £  are  unchanged,  althouah  the  convergence  to  the 
asymptotic  orbits  may  be  exceedingly  slow.  Hence  Theorem  £ 
remains  valid  if  we  allow  B  to  be  defective,  and  require  that 
V  have  a  nonzero  projection  on  a  principal  vector  C'f  B 
associated       with       an       eigenvalue       of        maximum       modulus.  If 

J=R  »*B*R  is  the  Jordan  normal  form  of  B  with  the  dominant 
eigenvalue  C  of  multiplicity  p!)  in  the  first  blockCsJ  along  the 
diagonal ,     then    these    vectors    have    the    form 

R*e 

I 

C  i  =1  , .  .  .  , pD ,     where    <e  >  is    the    natural    basi  s    f or    R    . 

5.     The    Power    Iteration 
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method  for  the  numerical  determinaiion  of  the  modulus  of  the 
dominant  eigenvalue  of  a  real  matrix  when  the  resulting  orbit 
is  periodic.  If  the  orbit  has  period  n>l  ,  then  the  modulus  is 
approximately  the  geometric  mean  of  the  Euclidean  norms  of  n 
consecutive  iterates.  Additionally,  the  argument  is  such  that 
the  nt_h  power  of  both  the  eigenvalue  and  it's  conjugate  are 
positive  real  numbers,  so  that  the  desired  eigenvalue  \  is 
gi  ven  by 

X  =  \X  |  w  C6.  13) 

for   som*-'   i  =1  ,  .  .  .  ,  n,   where   to  ,  .  .  .  ,  to   are   the   nth   root  s   of 

1  n 

unity.  If  the  orbit  is  aperiodic,  the  dominant  eigenvalues 
are  comp  ex  conjugates  and  fail  to  satisfy  C6.1D. 

As  an  example  of  the  use  of  the  power  iteration,   consider 
the  mat r i x 

which  hr-3  eigenvalues  ±V2.  In  114  p.  98]  it  is  shown  that  the 
powfr   method  CI.  21)   applied   to  this   iteration  Cwith  initial 

T 

vector  (  1,1!)  D  settles  into  a  period  two  cycle  on  the  two 
vectors 

i.e.  thfe  method  fails.  The  power  iteration  CI. ID  applied  to 
this  matrix  yields  asymptotically  the  period  two  cycle  Cfrom 
the    same    initial    vector) 

Since  this  is  an  asymmetric  period  two  cycle,  from  Theorem  2 
it  follows  that  there  are  both  positive  and  negative  dominant 
eigenvalues,     and    that    the    modulus    of    these    eigenvalues    is 


IX  I     =    C  II x   II  •  II x    II D1'"2 

1  2 

=  t  c  vs  /-/si  ■  c  z-/z/V5j  ] 1XZ 

Hence  the  eigenvalues  of  A  are  ±V2,  as  expected.  Similar 
results  would  be  obtained  using  any  norm  in  place  of  the 
Euclidean  norm  in  CI. ID. 

7.  Conclusions 

The  power  iteration  is  certainly  slower  than  the  power 
met  hoc  C  due  to  the  need  to  calculate  a  Euclidean  vector  no;  m 
rather  "  nan  simply  locating  an  element  of  the  vector  with 
maximum  modulus) ,  and  additional  information  can  be  gained  in 
only  a  restricted  set  of  cases  C primarily  when  the  dominant 
eigenvalues  are  complex  conjugates  and  real  multiples  of  a 
root  of  unity).  Nonetheless,  in  tnese  cases  it  does  provide 
useful  information  about  the  eigenvalues,  and  in  more  general 
case;  i  may  provide  some  insight  as  well.  For  example, 
inspection  of  the  elliptical  orbits  C Theorem  1)  can  be  used  to 
provide  information  about  the  existence  and  multiplicity  of 
null  eigenvalues.  For  these  reasons  the  power  iteration  may 
be  useful  in  certain  circumstances.  Of  course,  the  power 
iteraticn/power  method  is  in  some  sense  the  basis  of  most 
iterative  methods  for  the  eigenvalue  problem  [15]  and  so  this 
analysis  may  be  useful  in  the  analysis  of  more  practical 
algorithms  for  this  problem. 

We  mention  that  computer  gr  a  phi  ca.j  evidence  seems  to  indicate 
that  the  unsigned  power  method 


-     /  _^-.       .  i;  ,  .     i  r  t    -i  " 

y  ~      t->>*y     '    uy      n  ^    .    .   a  - 

r>-»l  r*  n     00 

has     the     proDerty     that      all      iterates      lie     on     a      /        conic-      C  wj  t  h 

oo 

respect  to  some  rotation  of  the  axes3  asymptotically.  As 
noted  in  [&,  p.  3621  ,  however,  the  use  of  the  Euclidean  norrr,  in 
CI. 13  greatly  facilitates  the  analysis  of  the  power  iteration 
Csee  also  [5,  p.  351  3  3  and  we  have  been  unable  to  show  a 
corresponding  result  for  the  unsigned  power  method  C7.13.  The 
analysis  of  cases  in  Theorem  £  did  not  depend  on  which  norm 
was  use  '  in  CI.  13  and  as  such  it  holds  for  C7.  ID  as  well  ; 
therefore  the  comments   in  §6,   excepting  those  about   the  t 

2 

conic  orbits,  a_e  equally  applicable  to  the  unsigned  power 
method,  which  is  no  more  costly  than  the  usual  power  method 
Crequinng  only  thai  an  extra  absolute  value  be  taken — bui 
t.  hi  s  would  be  done  during  the  corresponding  search  for  fj  3. 
Although  using  CT.13  in  place  of  C 1 . 23  means  that  the 
ai  goi  j  thi  will  not  converge  Cin  th  usual  sense"'  in  the  case 
of  a  do-  i  riant  negative  real  eigenvalue,  the  method,  properly 
inter  pre  ,-ed  Csee  also  [6,  p.  £573  3,  would  avoid  some  of  the 
problems     encouni  *--  e  i     when     using     the     power     method.  For      the 

matrix    A    given    in    C6    £3,     the    unsigned    power     method    C7. 13,     with 

T 

initial     vector    CI  ,13     ,     gives    a    period    two    cycle    on    the    vectors 

and    ag*in    the    fact    that    this    is    an    asymmetric    period    two    orbit 
implies     that    the    eigenvalues    of    A    are    ±  j  X.  |  ,     where 

|X|     =    C  II  x   II     •  II  x    II    31/Z 

1      00  2     00 

=    Cei  •  13 

=    V2 

as       expected.  This       type       of       reasoning       could       easily       be 
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